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Abstract 

We present a method to measure the magnetic properties of monolayers and ultra-thin films of magnetic material. 
The method is based on low energy muon spin rotation and y6-detected nuclear magnetic resonance measurements. A 
spin probe is used as a "proximal" magnetometer by implanting it in the substrate, just below the magnetic material. 
We calculate the expected magnetic field distribution sensed by the probe and discuss its temperature and implantation 
depth dependencies. This method is highly suitable for measuring the magnetic properties of monolayers of single 
molecule magnets, but can also be extended to ultra-thin magnetic films. 

Keywords: Single molecule magnets, ultra-thin magnetic films, monolayer, proximal magnetometry 



1. Introduction 

Recent developments of low energy muon spin rotation (LE-//SR) |[T1 |2l and /3-detected nuclear magnetic reso- 
nance (y6-NMR) |T/4l, provide unique local spin probe tools for measurements in thin films and multilayers. However, 
the application of these methods is limited by a minimal thickness of films, which provides sufficient stopping power 
for the implanted probes. For typical density materials a minimal 1-2 nm thickness is required at the lowest available 
implantation energy (1-2 keV). Therefore, these methods are generally not suitable for studies of monolayers and 
ultra-thin films. Nevertheless, for magnetic materials it is possible to perform a "proximal" measurement by implant- 
ing the probe in the substrate (or an under-layer) of the material, and sensing the dipolar magnetic fields from the layer 
of interest |5 , 6 |. 

To date, measurements of the magnetic properties of monolayers of single molecule magnets Q (SMMs) have 
been mostly performed using X-ray absorption spectroscopy (XAS) and X-ray magnetic circular dichorism (XMCD) 
imm. These are typically limited to high magnetic fields and provide information regarding the average static mag- 
netic properties of the irradiated portion of the monolayer In contrast, LE-yifSR and y6-NMR provide local probe 
measurements with important advantages; sensitivity to a large range of spin fluctuations and dynamics and applica- 
bility in zero or any applied magnetic field. However, a detailed interpretation of the measured spectra using LE-yuSR 
||6l or y6-NMR |5 1 when implanting the probe just below the monolayer is still absent. In this paper we calculate the 
magnetic field distribution in a non-magnetic substrate due to a monolayer of SMMs. This distribution can be used 
to fit LE-yuSR and /?-NMR measurements, providing information regarding the magnetic and geometric properties of 
the monolayer. The calculations can be easily generalized for the case of ultra-thin magnetic films by considering the 
domains as individual magnetic moments with a given average size. 
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2. Uniformly Magnetized Slieet Approximation 

As a simple approximation, a monolayer of magnetic moments (or a thin film) can be viewed as a uniformly 
magnetized sheet. Assuming magnetic moment M per unit area, located in the plane z - (see inset of Fig.[T]i with 
the magnetization aligned along -t The magnetic scalar potential from an annular region of radius p and width dp is 

Mzpdp 

d(l>M(p,z)^-^^-^—^. (1) 

The contribution of such a region to the magnetic field at a distance z from its center is 

dB,ip,z) = = -^ (^2^, 2)5/2 (2) 

where jUq is the permeability of the vacuum. We plot this contribution as a function of piz in Fig. [T] Integrating the 
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Figure 1: The contribution dB- of one annular region as a function Figure 2: Tiie z component of the field, sensed by a spin probe at 

of pIz- The inset shows a uniformly magnetized sheet divided into z = a/2, as a function of position. The color scale represents the 

annular regions. strength of the field. 

contributions from all annuli gives a magnetic field zero as expected. 

However, this approximation fails when looking at the microscopic structure of the sheet, i.e., as a collection of 
individual magnetic moments or even randomly oriented magnetic domains. This is the case when measuring the field 
due to such sheet using a local spin probe, which is positioned at a distance smaller than the microscopic details of the 
sheet (distance between moments or size of domains). For example, if we have a monolayer of moments with average 
spacing a between them, then the magnetic field sensed by a spin probe at a distance z < a from the monolayer is 
non-zero. 



3. Monolayer of Magnetic Moments 

In order to adequately account for the microscopic structure of a monolayer of magnetic moments, //,, one needs 
to sum the contribution of all moments in the monolayer. The dipolar field from fii at ri - XjX + yiy + ZiZ is 

Biird = ^ . (3) 

An rf 

A typical NMR or transverse field /iSR measurement is sensitive to the field distribution of one component of the field. 
In what follows we will restrict ourselves to the B, component, where z is normal to the monolayer surface. Working 
in polar coordinates with /?,■ - fi'^p + fxlz and ri - pip + ZiZ we can write 

p/, , //Q 3ziP■,^i'p + 2z^^^il-p',^il 

BXpi,Zi) = -J T--- . (4) 

47T (p2 + ^2)5/2 
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Therefore the field distribution, n(z,B), can be obtained by summing the contribution from all magnetic moments. 
For simplicity, we start by considering a monolayer of moments //, = fiz in the plane z = 0, arranged on a triangular 
lattice with lattice constant a. For example, taking /i = fig, a spin probe at z = a/2 will sense dipolar fields as shown 
in Fig.|2] The resulting field distribution is shown in Fig. [Sj^a), where different curves are results of summation of a 
different number of moments, such that is the radius (in units of a) of the considered "patch". In order to sample a 
field distribution as in Fig.|2j it is sufficient to sample only a unit cell of an equilateral (a) triangular area between 3 
moment^ Note that the qualitative shape of the distribution is almost independent of A^. However, due to the finite 
number of considered moments, we find a strong shift which decreases as a function of A^. 




Ba [J A'] 

Figure 3: The distribution of B, sensed by a spin probe due to a patch of moments with radius N (a) before and (b) after the correction in Eq. |6] 
The inset shows the most probable field as a function of N before (full symbols) and after (empty symbols) correction. 



In order to correct for this "artifact" we need to account for all the moments in the infinite monolayer However, 
since the moments outside the patch satisfy the criterion of a uniformly magnetized sheet, i.e. r,- » a, their contribution 
can be calculated as an infinite sheet of uniform magnetization with a hole of radius po - ciN around the origin. This 
can be calculated by integrating dB- (Eq. |2]) between p - pa and infinity. 



2p„(l+ ^2)3/2 



AB,= dB,(p,z)^ ^ r \o..n ^ (5) 



with ( = z/po- In contrast to the case of a uniformly magnetized sheet, AB, here does not vanish for finite po- Note 
that when the spin probe is far from the monolayer, such that po » z, then AB, « . This is quite different from 

the three dimensional case, where the remaining contribution is (where in this case M is the three-dimensional 
magnetization, i.e. the magnetic moment per unit volume), independent of how far we integrate out to. 

Assuming uniform magnetization due to moments fi = fig, ordered on a triangular lattice with lattice constant a, 
the magnetic moment per unit area is M = iibK V3a2/4). Using - po/a we arrive at, 

AB- . . (6) 

V3fl3A?(l +^2)3/2 



Due to the symmetry of the lattice the sampling can also be done on 1/6 of this area 
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This is the cause of shift seen in Fig. |3ja), and should be subtracted from the dipolar field of the patch of moments 
to correct for its finite size effect. Comparison between the field distribution before (Fig. |3ja)) and after (Fig. [3jb)) 
the correction shows clearly how well it works. Taking into consideration a patch with = 10 is already sufficient to 
accurately calculate the field distribution. Note that without the correction, even a patch with N = 60 exhibits a small 
artificial shift (inset of Fig. [3]). 

The field distribution has a similar shape to that in the vortex state of a superconductor IfTOll . The high field cut- 
off' is due to probes stopping exactly below a magnetic moment, while the low field cusp is due to probes stopping 
in middle between three neighbouring moments. However, one important difference between this case and a vortex 
lattice is that in the former the field distribution depends strongly on the stopping depth of the probe while in a 
superconductor the depth dependence is much weaker IfTTlfTSl . As we show below, this strong depth dependence will 
alter the shape of the field distribution entirely in a real measurement. Nevertheless, it is important to point out here 
that the shape of n(z, B) is determined by the lattice constant a and and the size of magnetic moment yU. Hence, these 
parameters can, in principle, be extracted directly from a measurement of the field distribution. Also, the details of 
the lattice, e.g. triangular or other, do not alter the qualitative shape of niz, B) dramatically. 



4. Measured Field Distribution 



In a LE-/vSR or y6-NMR measurements, the implanted probes stop in the sample at a distribution of depths, D{z), 
rather than a unique single depth. This distribution can be simulated accurately IIT3lfT4l using Trim.SP Monte-Carlo 
code ITSl . Therefore, the calculated n{z,B) should be averaged over D(z). Moreover, since the probe stops in the 
substrate of the monolayer, its intrinsic line shape (field distribution) has to be convoluted with the depth averaged 
field distribution. Finally, disorder in the monolayer, e.g. deviations from a perfect triangular lattice, can be taken 
into account by convoluting the depth averaged distribution with a Gaussian of width that represents the degree of 
disorder. In practice, one can simply convolute the calculated distribution with a single Gaussian or Lorentzian that 
includes both substrate and disorder effects. The resulting field distribution is then 

L(B)^I(B)* I D(z)n(z,B)dz, (7) 

•Jzo 

where 1(B) is the broadening due to the substrate and disorder, and zo is the distance between the magnetic cores in 
the monolayer and the surface of the substrate. In the case of SMMs, zq is due to the ligands used to graft the magnetic 
corse of a SMM to the substrate. In Fig. |4] we plot a typical field distribution due to a monolayer of Ifis moments, 
with a ~ 5 nm, zo ~ 1 nm and a stopping profile as shown in the inset. Both, the distribution with (solid line) and 




Figure 4: The field distribution L{B) sensed by a spin probe due to 
a monolayer of moments at z = 0. The dashed and solid lines are 
calculated with and without broadening, respectively. The inset 
shows a typical stopping profile of implanted probes in the sub- 
strate (shaded area). 



Figure 5: The square root of the second moment of the field dis- 
tribution as a function of distance from the monolayer (circles and 
solid line in inset). The inset is an expanded view of the shaded 
area. The dashed line is the fit described in the text. 



without (dashed line) the broadening 1(B) are shown. Surprisingly, the asymmetry of the distribution after averaging 
over depth is opposite to that of a single depth, i.e. we find a low field instead of the high field tail in Fig. [3] This is 
due to a larger near-zero field contribution of probes implanted deeper into the substrate, far away from the magnetic 
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cores in the monolayer. Note that unlike Fig. [3] the most probable field in L{B) is much closer to B = 0, and therefore 
in a real measurement we do not expect to see a considerable shift in the precession frequency of the probe, instead 
only some broadening can be measured [5ji6J. 

Now we turn to the depth dependence of the field distribution. As we mentioned earlier, we expect that for a depth 
z » fl the dipolar field of the magnetic moment vanishes as in the case of a uniformly magnetized infinite sheet. The 
depth dependence can be very useful in measurements as a tool for identifying (I) the nature of the sensed magnetic 
fields by the probe, dipolar or other, and (II) to estimate the characteristic length scale, a, in the monolayer [51. In 
Fig.|5]we plot the second moment of the field distribution as a function of depth. The dashed line is a linear fit of the 
second moment as a function of (a/z)^ in the range 50-100. Clearly, the second moment is inversely proportional to 
z^. This is in contrast to the 1 /z"^ dependence of the dipolar field from a single moment. A very small deviation from 
the (a/z)^ behaviour can be seen for z > a/3 (see inset of Fig.jsjl. Note also that the second moment vanishes at z ~ a, 
as expected. 

5. Information Drawn from an Experiment 

Considering actual measurements, in LE-juSR we measure the time dependence of the polarization of the im- 
planted muons in an applied transverese field. This is proportional to the Fourier transform of field distribution sensed 
by the muons, so that the average precession frequency is propotional to the average field and damping rate is propor- 
tional to the width of field distribution. In a j6-NMR experiment we measure the '^Li nuclear magnetic resonance line 
which directly proportional to the field distribution. In both cases it is, in principle, possible to fit the measured spec- 
tra. However, the asymmetry of the field distribution is quite small, and therefore hard to detect within the uncertainty 
of the measurement. Usually it is sufficient to approximately fit the data assuming a simple Lorentzian field distribu- 
tion. As we mentioned above, there is no detectable shift in the average field as a function of depth and temperature. 
In contrast, the width of the distribution depends strongly on these parameters. In fact, the width is proportional to 
il-P'}. Therefore, its temperature dependence can be obtained by measuring the temperature dependence of the field 
distribution width at a fixed depth (implantation energy). From this one should subtract the additional broadening due 
to the substrate. This can be easily evaluated by a measurement of the field distribution as a function of temperature 
dependence at large depths (high implantation energy). Finally, we point out that since we are dealing with local probe 
measurements, the averaging in (yU^), is on individual moments during the lifetime of the probe. This is different from 
the case of a XMCD measurement, where the averaging is done on the absorption signal from all molecules. 

6. Conclusions 

We have presented a detailed calculation of the properties of the magnetic field distribution, measured in the 
substrate, due to a monolayer near its surface. This distribution and its depth dependence can be used to infer the 
properties of the monolayer, such as the size of magnetic moment and the average distance between neighbouring 
moments. Although we restricted our calculations to oriented moments on a triangular lattice, they can be easily 
generalized for different systems. Nevertheless, the depth dependence of the width of the distribution (i.e. relaxation 
rate of the muon spin precession) and its temperature and depth dependencies should not be affected by the geometric 
details of the lattice. 
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